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This study evaluates the accuracy of an established 
reliability measurement procedure (NAVWEPS OD 29304) by 
computer simulation. The reliability measurement procedure 
assumes components fail according to an Exponential Failure 
Law. This study tests the accuracy of that procedure when 
components obey a m/eibull Failure Law or a Log Normal Fail- 
ure Law. 
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CHAPTER I 



INTRODUCTION 

The highly complex weapon systems being developed 
for the Navy today require a rapid assessment of system and 
subsystem reliability during research and development, and 
production phases. Currently there exists a statistical 
model, the Guide Manual For Reliability Measurement program 
(NAVWEPS OD 29304), that "...can be utilized by all contract- 
ors for subsystem reliability measurement and by the Navy 
for weapon system reliability measurement. This model is 
versatile in that it permits the combination of test data 
from all levels and can be continuously updated as new test 
data becomes available. The model is thus a rapid approx- 
imation procedure for determining system reliability. How- 
ever, the procedure has some restrictive assumptions which 
must be kept in mind. 

One notable restrictive assumption is that all compo- 
nents are assumed to have a constant failure rate. This may 
not be necessarily true, although many times is a good 
approximation. This leads to the commonly used Exponential 
Failure Law. However, suppose the failures were actually 
given by a Wei bull or a Log Normal Failure Law. Then this 
procedure must adequately estimate the system or subsystem 

■'■Guide Manual For Reliability Measurement Program 
(NAVWEPS OD 29304, 15 May 196^), p. 3-1. 



7 



reliability under the Exponential Law assumption. 

The purpose of this study is to evaluate the accuracy 

of the procedure in NAVNEPS OD 2930^ when failures actually 

occur in a i^eibull and Log Normal fashion. 

The approximation procedure uses statistical estimates 

of failure rates based on sample data and are thus subject 

to statistical uncertainty. Therefore, at best, the pro- ^ 

cedure must yield a lower confidence limit on system or 

subsystem reliability, or an upper limit on failure rate. 

It is desired to test the accuracy of this procedure by 

simulating the distribution of the lower confidence limit 

for various systems and comparing the simulation results 

with the true reliability. 

The method of simulation involves obtaining random 

Weibull and Log Normal failure time variates for components 

in a given system with a true system reliability, R . 

s 

These quantities represent the times to failure of compo- 
nents provided they are less than the duration of the test 
(i.e. the planned test time). Then, using the statistical 
procedure in the Guide Manual For Reliability Measurement 

Program, a 100(1 -&)% lower confidence limit, R_ T for 

s,na) 

R s is obtained. A large number of replications of this 
process are completed to construct a distribution of R g ) 
If the procedure is indeed accurate, the (l-a) t ^ 1 percentile 
point of the constructed distribution will be R c . That is, 

S 

A 

if ^s,L(a) a ^ rue 100(1 - &)% lower confidence limit for 
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If the (l-a) tlrl percentile 



R s> Chen p ( R s^ R S ,L(a)> = 1 - a - 

/\ 



point of the distribution of R Sj L( a ) denoted by A^_ a , 
then A]_ _ a should equal R g , regardless of the set of param- 
eter values chosen. Thus a measure of the accuracy of the 



procedure is given by the quantity 



A 1 - a “ 



Two other 



measures of the accuracy of this procedure are the mean and 

/\ 

standard deviation of the distribution of R S) L( a )» 

The accuracy of this procedure is examined for a 
variety of combinations of components and sets of parameter 
values. In all cases the systems under consideration con- 
sist of four components in series. All components are 
assumed to behave independently so the system reliability 
is the product of the four component reliabilities. 

As expected, the closer the deibull and Log Normal 
distributions approximate the Exponential, the better the 
accuracy of the procedure. However, for some components 
with failure time distributions differing considerably from 
the Exponential, the accuracy is fairly good using this 
method of reliability measurement. If the absolute differ- 
ence A-^ _ a - R S J is less than .02, the accuracy is called 
good. The accuracy of the procedure is good for planned 
test times of short duration and diminishes as they are 
increased. The study concludes that the reliability measure- 
ment procedure, being an approximation procedure, is useful 
for some systems following tfeibull and Log Normal Failure 
Laws provided they do not deviate from the Exponential Fail- 
ure Law by too significant an amount. 
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CHAPTER II 



STATISTICAL RELIABILITY MODEL 2 

The statistical Reliability Model provides the guide 
lines for reliability measurement on Navy weapons systems. 

The "Guide Manual" can be used by both contractors for sub- 
system reliability measurement and by the Navy for weapon 
system reliability measurement. A brief description of this 
model's assumptions and methodology are presented in this 
chapter. 

The system reliability obtained is based upon testing 
components for a time duration, called the planned test time, 
or until failure occurs under a given stress condition. 

Also, the model provides for testing a system for component 
failure by the same method. 

This model considers components in a complex system 
which must operate successfully over a defined "mission", 
lasting a specified duration of time. Time can thus be 
measured in "mission units". System reliability is thus a 
function of environmental and usage stresses which include 
vibration, shock, etc., and the operating vs. non-operating 
conditions as well as a function of the associated part 
failure rate parameters together with the time duration of 
the environments and usage stresses. The reliability 

^Guide Manual For Reliability Measurement Program 
(NAVWEPS OD 29304, 15 May 1965), PP. 3-1 to 3-15. 
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measurement system is dependent upon the following assump- 
tions : 

1. Constant failure rate. The exponential failure 

law is assumed to hold, mainly for its common 
usage, mathematical simplicity, practical 
simplicity, and reasonability. 

2 . Additivity of stress effects. "The failure rate 

induced by two simultaneously acting stresses 
is equal to the sum of the failure rates due to 
the two stresses acting sequentially. "3 Thus, 
mission experience can be simulated by adding 
data from separate environmental tests. 

3. Independence of component failures. "...This is 

assumed because components are normally tested 
individually by type, and subsystem reliability 
is estimated using component and other appli- 
cable test results. "4 

4. Failure rate constancy. "The failure rate is con- 

sidered a function of only the stress acting. "5 
In other words, the component being tested has 
no "memory" as to previous stresses of a differ- 
ent type. 



Computation of the system reliability is accomplished 
by the following procedure: 

An unbiased failure rate estimate is computed, pro- 
viding the failure rates are anticipated as being small. 
For ease of computation, a simplified unbiased estimator, 
which is slightly conservative, is used. For the i^ 11 com- 

A 

ponent, the unbiased estimator, A,^ , is given by: 



3 Ibid , p.3-3. 
4 Ibid . 

5 ibid. 
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( 1 ) 



A 




fi . 2 »1 

% 

2 T,, 2N 1 + 1 

j=l J 



th 

where = the sample size of the i component (i.e. the 
number of components of type i tested.),, 

T = the test time (time to failure) of the component 
of type i. 

T q ^ = the planned test time for component i, 
and f = the number of failure times which are less 
than T q1 . 

The system failure rate estimate is then given by: 



a m 

A A 

\ = 2 X, 

1=1 2 



( 2 ) 



where m is the number of components in the system. The 
component number assumes values i = l,„..,m. 

The variance of the unbiased failure rate estimate is 
given by: 



m 



, A 



Variance = 2 (1 /S ) 

i=l 1 1 



where 



N i 

S = 2 T. . = the sum of all test times 

i j=1 U 



accumulated on the N components of type i. 
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These estimates of failure rates for components In a 
given test condition can be continually updated during the 
development phase as more data becomes available. Also, 
during any higher level of assemblage, further testing can 
also yield data to modify the failure rate estimate. Data 
is obtained indicating the time that each component operates 
during the system test. 

The reliability estimate for the system is then given 

A 

A -X A 

by H = e where k is the appropriate failure rate esti- 
s 

mate and only series systems or subsystems are considered, 
given by formula (2). 

The reliability estimates obtained are subject to 
statistical uncertainty, therefore the important item to 
observe is the confidence interval about the estimate. 

Here, the lower confidence limit on system reliability 
becomes the pertinent quantity to observe. The "Guide 
Manual" bases the lower confidence limit on Normal theory 
and it is corrected to compensate for small values of X. 

The upper limit on failure rate, , provides the corres- 

u A 

ponding lower confidence limit on reliability by R. T = e~^ u . 
The upper limit on failure rate is given by: 




2a 

K c 




\ , a 2a 4a2 

V 4XK c + K c 
2 



K 



n 



m 

2 

1=1 



(1/S 1 ) 



m 

for 2 f . > 0 
1=1 1 

( 3 ) 

m 

for 2 f = 0 
1=1 1 
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where 



/\ 

c 



m 

£ 

i=l 



( V s i ) 

7 



n = the number of component-environment-test condition 
terms in the summation, and K = the percentage point of the 
Normal Distribution. The values of K for a given confidence 
level are not completely appropriate for small values of X. 
Consequently, a correction factor, 0, is used to obtain 
desired precision. Beta values are tabulated for the 80$ 
confidence limit in reference (2). Thus the percentage point 
of the Normal Distribution, K, is modified by 3K and replaced 
in formula (3). 

A more detailed description of this reliability meas- 
urement method plus examples of the procedure are given in 
reference (2). 
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SIMULATION PROCEDURE 

Suppose the assumption of constant failure rate used 
in the previously described statistical model is relaxed. 

It is now desired to test the accuracy of that model when 
failures occur according to a Neibull Failure Law and a Log 
Normal Failure Law. For testing purposes by simulation 
techniques, consider a series system with four components, 
each component having a well defined failure law. 

The true reliability for each component of this series 
system is given by R^ , i = 1 , 2 , 3 , 4, and thus the true system 
reliability is 

4 

R = 2 R . (4) 

s i=l 1 

th. 

Using the statistical model, an estimate of the (1-a) 

/\ 

percentile confidence limit of the system, R ,, can be 

s,na ; 

obtained by computer simulation. This is a random variable 
as calculated from the statistical model. If 

P( R g > R S} L( a ) ) = 1 “ a 

A 4-U 

then, in fact, R g is an exact (1-a) percentile lower 

confidence limit for R s . This says that R g is always the 
(l-a)^ percentile point of the probability distribution of 

A /\ 

R g L( a )* The distribution of R g is constructed by 
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“h V"» 

computer simulation. Letting be the (1-ct) per- 

m's 

centile point of the distribution of R g L( a )> then the 

absolute difference, |A-^_ a - R g J , is a measure of the 

accuracy of the statistical model. 

Two other measures of the accuracy of this model are 

the estimated mean, R g anc * the estimated standard 

deviation, s£ , which give some assurance that the 

K s,L(a) 

actual values generated by the procedure are reasonable. 

The simulation procedure itself is now discussed for 
the series system of four components: 

The distribution of l* s L( a ) constructed by gener- 
ating 500 random observations on f° r a Si yen system, 

number of components tested, and reliability of the com- 
ponents. For this study a four component series system is 
used, with each component having reliability R^ , i = 1,2,3,^. 
The reliability, R^ , is defined for each component by 
R^ = R^l) = P( T^ >1), where T^ is the time to failure 
random variable for a component of type i. This establishes 
the parameter values for a particular failure rate distri- 
bution. For example, if a component is to obey a Weibull 
Failure Law, 

R(t) = e _(u ) 3 (5) 

_ x 3 

and R(l) = e = . 995 j then, for 3=2, this implies that 
X = .0707. 

Similarly, for the Log Normal Failure Law, if 
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R(t) = e Z = 1 - § 



( 6 ) 



In t -/a 
cr 

and R(l) = 1 - (|>(-/'V<r ) = .995» then, for cT= 1 , this 
implies that /a = 2.576. The failure rate functions deter- 
mined from the reliabilities and selected parameters are 
displayed in Appendix I. 

A random number generator is used to obtain a sample 
of tested components of type i which fail according to 
their respective failure rate distributions. This Monte 
Carlo method draws a random number which is Uniformly Dis- 
tributed ( 0 , 1 ) and converts it to a Weibull variate or a 

Log Normal variate to obtain the time to failure of the i tn 

_ f \ m \ 3 

component, T . For the Weibull Distribution, e ' = Y, 

where Y is distributed as Uniform (0,1), implies that 



T = (-In t ) 17,3 
1 



(7) 



is a random Weibull variate. For the Log Normal Distri- 
bution, if X is a random variable distributed as Normal (0,1), 
then Z = crX + y*. is a random variable distributed as Normal 
(yA,<r^) and thus T = e z is a random Log Normal variate. In 
this manner, N^ failure times are generated for each of the 
four components ; 



T 11 T 12 


T 13 ° 


° • T 1 N; 


T 21 T 22 


T 2 3 * 


* • T 2N. 


t 31 t 32 


T 33 * 


* * T 3n 


T 4l t 42 


T 43 . 


• * T 4n 
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Using these failure times in the statistical procedure 

A A 

discussed in the previous chapter, values of R and R_ T , \ 
are obtained with that model. Each is examined for 

failure before the end of the planned test time or else 
terminated at the end of the planned test time. This same 
procedure is completed 500 times to generate 500 random 

A 

observations, thus constructing the distributions of R s 
and ®s,L(a)* 

The mean and standard deviation of the distributions 
are determined by: 



A _ ! 500 A 

R s,L(a) 500 R s,L(a) k 



( 8 ) 



a i 500 A 
R s - 500 iii R s k 



( 9 ) 



S R 

R s , L ( a ) 




(R s,L(a) k _ R s,L(a) ) 



( 10 ) 





500 

2 

k=l 




a . 2 

V 



(id 



The 500 values of ^( a ) are then ordered and the 
(l-a) tR percentile point of the distribution is obtained. 
Thus, the value is available for comparison with R_. 

This process is done merely by counting down to the 101 st 
value of R s from the largest, which represents the 

80^ percentile point. 
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Several cases are examined with various failure rate 
distributions and with various combinations of components. 
For each case examined, the planned test time, T Q ^, varies 
to produce an optimistic or pessimistic reliability test. 

If T 0 ^ = .5, the system will yield an optimistic test since 
the planned test time will be reached before most components 
fail. The opposite will occur when T Q ^ = 5. A T Q ^ chosen 
such that the average failure rate is equal to the failure 
rate at time one appears to be a good test of system relia- 
bility for this procedure. That is, make T Qi such that it 
satisfies the equation 




( 12 ) 



where z(t) is the failure rate function; the failure distri- 
bution divided by the reliability function: z ( t ) = f ( t )/R(t ) . 

All cases studied by simulation procedures are given in 
Table I. All simulation results of these cases are given in 
Table II. 
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TABLE I 



CASES STUDIED BY SIMULATION PROCEDURES 



Case 

No. 


Compo- 
nent No. 


Failure 

Law 


Parameters 


R i 


R 

s 


1 


1 


Wei bull 


3=1.5, a.= . 0293 


.995 


.98 




2 


Weibull 


3=2.0, A=. 0707 


.995 






3 


Log Normal 


cr=i. 0,^=2. 576 


.995 






4 


Log Normal 


CT=2.0,yM=5.152 


ON 

ON 

• 




2 


1 


Weibull 


3=1.5, >,= .0293 


.995 


.98 




2 


Weibull 


0=1.5, A=. 0293 


.995 






3 


Log Normal 


<r=2.o,/4=5.i52 


.995 






4 


Log Normal 


Q-=2.0,ylA=5.152 


.995 




3 


1 


Weibull 


3=1.33, 1=. 0378 


O- 

OO 

ON 

• 


.95 




2 


Wei bull 


0=1.33, ^=.0378 


.987 






3 


Log Normal 


cr =1.5, /^=3.345 


.987 






4 


Log Normal 


cr =1. 5,^=3. 345 


O- 

00 

ON 

• 




4 


1 


Weibull 


3=1.2, i=.0485 


.97^ 


.90 




2 


Wei bull 


3=1.2, i=. 0485 


.974 






3 


Log Normal 


or=1.5, y u.=2.919 


.974 






4 


Log Normal 


<1=1. 5, =2. 919 


.97^ 




5 


1 


Weibull 


3=1.2, a=. 0485 


.974 


.90 




2 


Wei bull 


3=1.2, \=.0485 


.974 






3 


Wei bull 


3=1.2, A.= . 0485 


.974 






4 


Weibull 


0=1.2, A=. 0485 


.974 
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TABLE I (Continued) 



Case 

No. 


Compo- 
nent No. 


Failure 

Law 


Parameters 


R i 


R s 


6 


1 


Log Normal 


<r=i .5, >*=2.919 


.974 


.90 




2 


Log Normal 


cr =1 . 5 » >*=2.919 


.974 






3 


Log Normal 


0"=l • 5 » >*=2 .919 


.974 






4 


Log Normal 


cr=l.5*>*=2.9l9 


.974 




7 


1 


Weibull 


3=1.1, A=. 0483 


.965 


.867 




2 


Wei bull 


3=1.1, *=.0483 


.965 






3 


Weibull 


3=1.1, A=. 0483 


.965 






4 


Weibull 


3=1.1, A=. 0483 


.965 




8 


1 


Log Normal 


C=2 . 0 , >*=4.0 


.977 


.912 




2 


Log Normal 


(T= 2.0, >*=4.0 


.977 






3 


Log Normal 


O' =2 . 0 , >*=4.0 


.977 






4 


Log Normal 


<T= 2.0, >*=4.0 


.977 
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CHAPTER IV 



RESULTS AND CONCLUSIONS 

In the preceding chapters the statistical procedure 
from the Guide Manual For Reliability Measurement Program 
has been briefly explained to obtain a 100 ( 1 - a )% lower con- 
fidence limit on system reliability. The simulation of this 
method on a computer has been explained and the technique 
for measuring the accuracy of the statistical procedure 
stated. Now, from the results of the simulation, some con- 
clusions may be drawn. 

The accuracy of the statistical procedure has been 
examined for eight cases which represent different combi- 
nations of components for various sets of parameter values. 
The results for these cases are given in Table II. The 
mean and standard deviation of the 500 values of R g ^( a ) 
and S s are given, along with the (l-a) tlrl percentile point, 
A]_ _ a . In all cases studied, a was taken to be .20. 

The accuracy of the procedure is measured by the 

t Vi 

absolute difference between the 80 percentile of the dis- 

A 

tribution of R g ^ 20)’ A 0 q, and the actual system 

reliability, R_ (i.e. accuracy = A Q n - R ). 

s I . oO s I 

Also stated in the results, is the quantity TT, which 
is the amount of testing relative to the component un- 
reliabilities. In other words, this is a measure of the 
amount of testing required to achieve a desired accuracy 
for a given system with component failure rates, z^(l), 
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and planned test times, T Qi . Then TT is given by 



m 

TT = 2 N. z. (1 )T . . 

i=l 11 01 



(13) 



As an example of the accuracy study, consider Case 1, 
where the system consists of four components in series dis- 
playing the following failure laws: 



Component 1 
Component 2 
Component 3 
Component 4 



Neibull, (3 = 1. 5» \=.0293 
•ieibull, (3=2.0, \=.0707 
Log Normal, <T= 1.0, 2.576 

Log Normal, <r = 2.0, jx= 5.152 



The component reliabilities, R^ = .995, for all four 

components, so the system reliability, R = . 98 OI. The 

s 



T 01 = .5, IS 



A .80 “ R s 



measure of accuracy for the optimistic planned test time, 

= .008 for 500 components of each 
type tested. The pessimistic planned test time, T Q ^ = 5.0, 
yields an accuracy of .055 for 100 components of each type 
tested. The ad hoc planned test times, T Q ^ = 2.25, T 0 £ = 2.0, 
Tq-j = 1.95, and T q ^ = 3.25, as determined from formula (12), 
give an accuracy of .030 for N^ = 100. The accuracy for 
the ad hoc planned test times is only fair. This planned 
test time represents the best measure of accuracy of the 
three selected. 

It can be seen from Table II and Appendix I that the 
accuracy of the procedure is a function of how closely the 
Weibull and Log Normal failure rate functions approximate 
a constant failure rate function. The greater the Weibull 
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and Log Normal failure rate functions deviate from the 
constant failure rate, the more inaccurate the procedure. 
Also, the accuracy decreases as the planned test time is 
increased. The Weibull failure rate function increases in 
time and the Log Normal failure rate function increases 
then decreases in time, thus deviating from the constant 
failure rate for large planned test times. 

It can thus be concluded from the study that the 
statistical procedure can be useful for reliability approx- 
imations of a system, if the components that fail according 
to these non-constant rates do not deviate too significantly 
from a constant failure rate. 
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SIMULATION RESULTS 
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APPENDIX I 



GRAPHS OF COMPONENT FAILURE RATE FUNCTIONS 
FOR EACH SYSTEM TESTED 




FIGURE 1 

Component Failure Rate Functions, z(t), for Case 1. 
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FIGURE 2 

Component Failure Rate Functions, z(t), for Case 2. 
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FIGURE 3 



Component Failure Rate Functions , z(t ) , for Case 3 * 
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FIGURE 4 

Component Failure Rate Functions, z(t), for Case 4. 
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FIGURE 5 

Component Failure Rate Functions, z(t), for Case 5* 
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FIGURE 6 

Component Failure Rate Functions, z(t), for Case 6* 
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FIGURE 7 

Component Failure Rate Functions, z(t), for Case 7. 
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FIGURE 8 

Component Failure Rate Functions, z(t), for Case 8. 
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on 



AP® END I X II 



COMPUTER PROGRAM 



01 MENS I ON PS (500) , RSL ( 50C ) . S ( 4 ) , NF ( 4 ) , 

♦ N( 4 ) ,HLAM(4) f BTA(30),TO(4 t 3) 

C RF AD IN SAMPLE SIZE, N 

READ 1 ,N 

1 FORMAT (418) 

C READ IN BFTA CORRECTION FACTORS 

READ 2, RTA 

2 FORMAT (5FR.3) 

C READ IN PLANNED TEST TIMES 

R c AD 3 » ( (T0( l f j) , i=i ,4) , J = l f 3) 

3 FORMAT (4F10.0) 

C R F AD IN FAILURE DISTRIBUTION PARAMETERS 

READ 12, ALFA1 ,BETA1, ALFA2, RFTA2, 

♦SIGMA! ,XMFANI , SIGMA2, XMEAN? 

1 2 FORMAT ( 8F6.0) 

PRINT 4 

4 BORMAT ( 1 H 1 ) 

C RUN SIMULATION FOR THREE DIFFEPENT PLANNED 

C TEST TImps 

DO 999 I RUN=1 , 3 

C URN IS RANDOM NUMBER GENFRATOR 

Y = (JRN ( 0 ) 

RTOT=C. 

RLTOT=D, 

C GENERATE 500 VALUFS 0*= LOWER CONFIDENCE LIMIT 

DO 50 K=1 ,BO n 
hlamd=o. 

FT0T=0. 

TS = 0. 

TSUM=0. 

C SIMULATE FAILURE OF FACH OF THF FOUR COMPONENTS 

DO 3 n 1=1,4 
S( I )=o. 

NF ( I ) =0 

M=N ( I ) 

C GENERATE N FAILURE TIMES FOR FACH COMPONFNT 

DO 20 J=1,M 
Go TO (5 ,4,7,8) , I 
c Y=URN( 1 ) 

B T NV1 = 1 • /BET A1 

T P= ( (-ALOG(Y) ) **R I NV1 ) / AL F A l 
GO TO 1C 

6 Y = U R N ( J. ) 

BTNV2=1. /RETA2 

TP = ( (-ALOG(Y) ) **BI NVR ) / AL F A2 
GO TO 10 

7 TSi=o. 

DO R I J=1 , 1 2 
<= T S 1 = T SI ♦•JR N ( 1 ) 

X! =TSl-6, 

71=SIGMAI*X1+XMEAN1 
TP=EXP( 71 ) 

GO TO 10 

8 TS2=0, 

DO 11 I J=i , 1 2 
11 TS2=TS2+URN(! > 

X2=TS2-6, 

Z2=SIGMA2*X2+XMEAN2 
TP=FXP( 72) 

10 CONTINUE 

COUNT NUMBER OF FAILURES BEFOPF PLANNED 
TEST TIMP EXPIRES 
IPfTP.LT.TOH, IRON) > GO TO 15 
T=TO( I , IRON) 
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on on o on 



GO TD 18 
15 T = TP 

NF( I>=NF<I > + l 

18 sm = sc i h- t 

20 CONTINUE 

C C4LCULATF LAMOA HAT ANO $ I MUL4 T ED SYSTFM 

C RELIABILITY 

XN = N < I ) 

F = NF( I ) 

HLAM(T)=(F/S(m*( 2.+XN/12. *XN+1. > ) 

FTOT=FTOT+F 
ML A MO =HL AMD+HL AM ( I ) 

tsum =tsum+(hlamcu/su n 
T$=TS+1» /SC I ) 

^0 CONTINUE 

RS< K1=EXP(-Hl AMO) 

RTQT =RT0T*RS(K) 

CALCULATE UPPER LAMOA HAT AND 
LOWER CONE I DFNCF LIMIT RFLIARILITY 
iftot=ftot+ 1 . 

TFC IFTOT-3! ) 31 ,32,3? 

M 8ETA=BTA( I FTOTl 
00 TO 33 

32 RFTA=1. 

33 x* = ( . R42»*RETA 
IP(FTOT) 35,35,4^ 

35 ULAMO=( X***2 )*TS/4. 

GO Tn 41 

4C chat=tsum/hlamd 

IJL A M0= ( 2 . ♦HLAMO+ ( XK**2) *CHAT + SQRT < 4. *HL AMO* 
*(XK**?>*CHAT+(XK**4>MCHAT**2» )>/?, 

M RSL(K»=EXP(-ULAMD) 

RLTOT=RLTOT+RFL(K ) 

50 continue 

CALCULATE MEANS ANO STANDARD DEVIATIONS 
RR AR = RTOT/^C r », 

S'|M »0 . 

DO 60 K = 1 » c 00 
0IFF=(RS(K»-PPAR|+*2 
SHM=SUM+DIFF 
6 C CONTINUE 

SR = SO°T ( SUM/500. ) 

RLRAP=RLT0T/5 nr '» 

$ijM=n. 

00 61 *=1 , 50^ 

OIFF=(RSL( K1-RLB4R>**2 

SUM = S(JM + OIFF 

61 CONTINUE 

SRL=SQRT(SUM/50O . ) 

OPOFR 500 VALUES OF LOWER CONFIDENCE LIMIT 
ON RELIABILITY ANO S C LECT ROTH PERCENTILE, A 
11 = 1 

00 80 K= 1 , 1 01 
B=RSL(l> 

00 70 1=1,500 
IF ( B-PSL ( I 1 ) 65,70,70 

65 R=RSL(I> 

1 I = I 

70 CONTINUE 
RSL( 111=0. 

4 = 8 

pr CONTINUE 

PRINT OUT VALUES OF MEANS AND STANDARD DEVIATIONS 
AND RTTH PERCENTILE 
PRINT 90, RLRAR ,5RL,RBAR,SR,A 
or criRMATl //1CX,5F2% 5) 

Poe CONTINUE 
ENO 
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APPENDIX III 



WEI BULL DISTRIBUTION 6 

The v/eibull Probability Distribution is defined by: 



3 

f (t ) = 3 X | 3 t |3 ' 1 e" 0 <t<°° 

where X is the ’’scale" parameter, and 3 is the "shape" 
parameter. 

The Cummulative Distribution Function (CDF) is then 




and the failure rate function is 



z(t) = f (t)/R(t) = X^t * 3 ” 1 

where R(t) = 1 - F(t) = the reliability function. The 
Weibull Probability Distribution reduces to the Exponential 
Probability Distribution when 3=1. Increasing failure 
rates are obtained only when 3 > 1 . 

^Daniel R. Cox, Renewal Theory (London: Methuen and 
Co. Ltd., 1962) p. 21. 



44 



APPENDIX IV 



LOG NORMAL DISTRIBUTION 7 



The Log Normal Probability Distribution is given by 

1 

f(t) = — - — e ■" J 0 < t < oo 

YItt frt 

= e-^ 



p 

where X is distributed as Normal { y U,C' ). For computational 
purposes, it is beneficial to derive the Cummulative Distri- 
bution Function as 



F(t) = P( t < T ) 

= P( t<e X ) 

= P(ln t < X ) 

then fit) = F' (t ) = -L- 4)[ ln 

The failure rate function is given by: 



z(t ) 



f(t)/R(t) = 



1 - 



7 Paul L. Meyer, Introductory probability and Sta- 
tistical Applications (Reading, Mass: Addison-Nesley , 19&5) 
p. 18?. 



45 



INITIAL DISTRIBUTION LIST 



No. Copies 



1. Defense Documentation Center 20 

Cameron Station 

Alexandria, Virginia 22314 

2. Library 2 

Naval Postgraduate School 

Monterey, California 93940 

3. Chief of Naval Operations (OP 9^) 1 

Department of the Navy 

Washington, D.C. 20350 

4. professor W. Max Woods 10 

Department of Operations Analysis 

Naval Postgraduate School 
Monterey, California 93940 

5. LT Kenneth Alan Huffman, USN 1 

2710 Petty Road 

Muncie, Indiana 47304 



46 



UNCLASSIFIED 



Si’i'iinlx Cl a s si t u'.if ton 



DOCUMENT CONTROL DATA - R & D 1 


i o 1 ' a>in a ini. a l t i \ i t > (i orpomte author) 

Naval Postgraduate School 
Monterey, California 93940 


20. REPORT SECURITY CLASSIFICATION 

UNCLASSIFIED | 


2h. GROUP 1 


EVALUATION OF TIIE ACCURACY OF A 

RELIABILITY MEASUREMENT PROCEDURE 1 

USING SIMULATION TEC INI iUES 



I’aesis 



? AUTHORiSi (Fuxf name, middle initial , last name) 



r.enneth Alan Huffman, Lieutenant, United States Navy 



6 REPORT DATE 

December 1968 


7a. TOTAL NO. OF PAGES 

46 


7b. NO. OF REFS 

4 


CONTRACT OR GRANT NO 

b. PROJEC T NO 

c. 

d. 


9a. ORIGINATOR'S REPORT NUMB 


E R (S) 


9b ' thl s* report) P ° R T NO(S) ( An Y other numbers that may be assigned 


K DISTRIBUTION STATEMENT - ■ 

Distribution of this document is unlimited. I 


" ' SUPPLEMENTARY NOTES 
3 ABSTRACT 


12. SPONSORING Ml LI T ARY ACTIVITY 1 

Naval Postgraduate School 
Monterey, California 93940 



fnis study evaluates tne accuracy of an established 
reliability measurement procedure (NAVNEPS CD 29304) by computer 
simulation. The reliability measurement procedure assumes com- 
ponents fail accurding to an Exponential Failure Law. This study 
tests tie accuracy of that procedure when components obey a 
.veibull Failure Law or a Log Normal Failure Law. 
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